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Abstract 

The presence of non-Gaussian features in the Cosmic Microwave Background (CMB) 
radiation maps represents one of the most long-awaited clues in the search for the actual 
structure of the primordial radiation, still needing confirmation. These features could shed 
some light on the non trivial task of distinguishing the real source of the primeval pertur- 
bations leading to large scale structure. One of the simplest non-Gaussian signals to search 
is the (dimensionless) skewness S. Explicit computations for S are presented in the frame 
of physically motivated inflationary models (natural, intermediate and polynomial potential 
inflation) in the hope of finding values in agreement with estimated quantities from large 
angle scale (e.g., COBE DMR) maps. In all the cases considered the non-Gaussian effects 
turn out to lie below the level of theoretical uncertainty (cosmic variance). The possibility 
of unveiling the signal for S with multiple-field models is also discussed. 
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1 Introduction 



The existence of anisotropies in the CMB radiation as recently detected by CO BE []T]] and sub- 
sequently confirmed both by balloon-borne scans at shorter wavelength |^ and by ground- 
based intermediate angular scale observations ^ has triggered a large body of literature 
dealing with the non trivial task of finding the correct statistics able to disentangle the rele- 
vant information out of the primeval radiation maps. Both small and large angle scale probes 
of the microwave sky have been and are the real data our theoretical models must reproduce 
before we might call them viable. Discrimination between the two main theories for the 
origin of primordial perturbations, namely, whether these are due to topological defects |Q 
produced during a GUT phase transition or to early infiationary quantum fiuctuations 0, 
has by now become a difficult matter. 

In this paper we will work in the frame of infiation and, in particular, we will be mainly 
concerned with single-field infiationary potentials. Some of the most popular models are 
characterized by their simplicity and universality (such as quadratic and quartic chaotic 
potentials), by their being exact solutions of the equations of motion for the infiaton (like 
power-law and intermediate infiation) or by their particle physics motivation (as natural 
infiation with an axion-like potential). In contrast with these simpler models where one has 
just one relevant parameter more general potentials, with more freedom, were also considered 
in the literature. One example of this is the polynomial potential which for an adequate 
choice of the parameters was found to lead to broken scale invariant spectra on a wide range 
of scales with interesting consequences for large scale structure. 

One should also worry about initial conditions . While for single-field models the only 
effect of kinetic terms consists in slightly changing the initial value of (leaving invariant 
the phase space of initial field values leading to sufficient infiation), for models with more 
than one scalar field initial conditions can be important, e.g., for double infiation ||^ (with 
two stages of infiation each one dominated by a different inflaton field) leading to primordial 
non-Gaussian perturbations on cosmological interesting scales. Within the latter models, 
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however, the question of how probable it is that a certain initial configuration will be realized 
in our neighbouring universe should be addressed. More recently other examples of inter- 
esting multiple-field models with broken scale invariance have also been considered (see e.g. 
Ref. 0). Here all scalar fields contribute to the energy density and non-Gaussian features 
are produced when the scalar fields pass over the interfaces of continuity of the potential. 
Extension of the single-field stochastic approach developed in Ref. [1^ for calculating the 



CMB angular bispectrum generated through Sachs- Wolfe effect from primordial curvature 



perturbations in the inflaton in order to include many scalar fields is therefore needed |]TT 

The plan of the paper is as follows. In Section 2 a brief overview of some general results 
is given while in Section 3 we concentrate on trying to extract numerical values for the non- 
Gaussian signal (the dimensionless skewness in this case) predicted in the frame of three 
different inflationary models. Section 4 contains some general conclusions. 

2 The CMB skewness 

We will here briefly summarize the steps that lead to the calculation of the mean two- and 
three-point functions of the temperature anisotropies and in particular to the zero-lag limit 
of the latter, namely the skewness. 

As usual we expand the temperature fluctuation in spherical harmonics ^ (x; 7) = 
Z^^^i Yl,m=-i (^T{^)^i^r{l)i where x specifles the position of the observer and the unit vec- 
tor 7 points in a given direction from x. We represents the window function of the speciflc 
experiment. Setting Wo = Wi = automatically accounts for both monopole and dipole 
subtraction; for i > 2 one can take Wi ^ exp l)cr^ , where a is the dispersion of 

the antenna-beam proflle, which measures the angular response of the detector (e.g. ||12|| ). 
In some cases the quadrupole term is also subtracted from the maps (e.g. Jll); in this case 
we also set W2 = 0. The multipole coefficients aj^ are here considered as zero-mean non- 
Gaussian random variables whose statistics derives from that of the gravitational potential 
through the Sachs- Wolfe relation ^(x;7) = |$(a; + ro7), where ro = 2/Ho is the horizon 
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distance and i^o the Hubble constant. 

In the frame of the inflationary model, the calculations reported in Ref . [|T0| lead to general 
expressions for the mean two- and three-point functions of the primordial gravitational 
potential, namely ($(ro7i)$(ro72)) = (9QV207r) E£>o(2^ + l)^^(7i ■ 72)^ and 

($(ro7i)$(ro72)$(ro73)) = ^^-^s E (2j + 1)(2^+ 

X [P,-(7i ■ 73)P^(7i ■ 72) + ^j(72 ■ 7i)i^£(72 ■ 73) + ^i(73 ■ 7i)^£(73 ■ 72)] , (1) 

with a Legendre polynomial and where $3 is a model-dependent coefficient. These ex- 
pectation values are a statistical average over the ensemble of possible observers and can 
only depend upon the needed number of angular separations. 

The ^-dependent coefficients are defined by {Qf) = (^^+^) Q^C^, with Q = (Ql)^^"^ the rms 
quadrupole, and are related to the gravitational potential power-spectrum Pq>{k) through 
Ci = /o°° dkk"^ P<^{k)iJ{krQ) / /q°° dkk'^Pq,{k)i2{krQ), where is the £-th order spherical Bessel 
function. The rms quadrupole is simply related to the quantity Qrms~ps defined in Ref.[jl], |13|: 
Q = \^A7rQrms-PS /To, with mean temperature Tq = 2.726 ± O.OIK |T^. For the scales of 
interest we can make the approximation P^{k) oc A;""^, where n corresponds to the pri- 
mordial index of density fiuctuations (e.g. n = 1 is the Zel'dovich, scale-invariant case), in 
which case P, ^ we have = T + f - i) T (f - f ) / (r + | - f ) T (| + f )) . The 



equations above allow us to compute the angular spectrum, {oil^a^^*) = 5£i£2'^mim2 
and the angular bispectrum, 

where the coefficients 'H1\'2'^' = / rffi^F^';'!* (7)^^7(7)^^7(7) are only non-zero if the in- 
dices £i, rrii {i = 1, 2, 3) fulfill the relations: \ij — ik\ ^ ^ \^j + ^k\, ^1 + ^2 + ^3 = even and 
mi = 7712 + ^3- From the last equation we obtain the general form of the mean three-point 
correlation function for the temperature perturbations whose explicit form we will not need 
in the present work. Some simplifications occur for the CMB mean skewness {Cs{0)) for 
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which we obtain 



25(47r) 



-$3 E (2^1 + 1)(2^2 + 1)(24 + + CeA^ + CeAi] 



(3) 



where the coefficients J-'ij^e^is = (^tt)^^/ dfi^^ / dQ^Pi-^^j ■ ^')Pi^{j ■ ^')Pf.^{'^ ■ 7') may be 
suitably expressed in terms of products of factorials of £1, £2 and £3, using standard relations 
for Clebsch-Gordan coefficients. 



By working in the frame of the stochastic approach to inflation |17, 18 1 we are able to 



compute the three-point function for the inflaton field perturbation 5(f) [|T0|- Among the 
primordial perturbation scales we find those that stretched to super-horizon sizes approxi- 
mately 60 e-foldings before the end of inflationary era to re-enter the Hubble radius next, 
during matter domination, as energy-density perturbations (and thus affecting also the grav- 
itational potential $). During decoupling these scales are still greater than the horizon and 
therefore no microphysics alters the primevally imprinted signal they carry. The absence of 
non-linear evolution for these large angle scales makes physics simple and therefore highly 
predictive. 

The stochastic analysis naturally takes into account all the multiplicative effects in the 
inflaton dynamics that are responsible for the non-Gaussian features. Extra-contributions 
to the three-point function of the gravitational potential also arise as a consequence of the 



non-linear relation between $ and 5(j) When all primordial second-order effects are 

taken into account we get 





r(3-n)r(| + f) 




r (2 - 1) 





for the rms quadrupole, while for the "dimensionless" skewness S 
find m] 



(4) 



we 



32772 



An) 



(5) 



where we denoted with Xqq the value of the steepness of the potential X(0) = mpV'{(f>) /V{(f)) 
evaluated at 06o (the value of the inflaton 60 e-foldings before the end of inflation) and where 
I{n) is a spectral index-dependent geometrical factor of order unity for interesting values of 
n. 

It is also of interest here to calculate the accurate form of the spectrum of primordial per- 
turbations (e.g. by finding the value of the spectral index) at the moment the scales relevant 
for our study left the Hubble radius. Our expression S for the dimensionless skewness is ac- 
curate up to second order in perturbation theory (the first non vanishing order) and therefore 
we should calculate n at least to the same order. We will borrow the notation for the slow- 
roll expansion parameters from Ref . ||2l| . These are defined as e(0) = {H' {(fy / H {(fyy /At:, 
r^{(j)) = mlH"{(j))/{ATTH{(f))) and ^(0) = mlH'" {({))/ {An H '{({))) 0. In the slow-roll approxi- 



mation e and t] are less than one. The same is not true in general for ^ and this may cause 
consistency problems when this term is incorrectly neglected (see the discussion in [^). We 
will see below examples of this. 

Let Q| {Qt) be the contribution of the scalar (tensor) perturbation to the variance 
of the quadrupole temperature anisotropy. The complete second-order expressions for the 
tensor to scalar ratio and for the spectral index are given respectively hy R = Q^IQ% ~ 
14e [1 - 2C{r] - e)] and 



n = 1 - 2e 



2-2+4(C + l)e-(5C + 3)r7 + Ce(l + 2(C+l)e-Cr7) , (6) 



evaluated at ^ 06o- In this equation C = — 2 + ln2-|-7^ —0.7296 and 7 = 0.577 is the 



Euler-Mascheroni constant [231. 



It was realized |]2^, ^ that a positive detection of a non-zero three-point function 
(an in particular a non-zero skewness) in the temperature fluctuations on the microwave sky 
does not imply an intrinsically non-Gaussian underlying field. This problem is related to the 
so-called "cosmic variance" . Relevant for large angular scale fluctuations, this variance is in 
fact the strongest source of (theoretical) noise we have to deal with and essentially reflects 
the impossibility of making observations in more than one universe. One way to quantify this 
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effect is through the rms skewness of a Gaussian field (Cf (0))^^^^^, which we may express 
as {C^{0))Gauss ~ 3 j\dcosa{C2{Q:))^ . A convenient quantity to compare with S is the 
normahzed rms skewness (C*! (0))^^f„^g/ {C2{0))^^'^- For interesting values of n this ratio is of 
the order ~ 0.1. A rough criterion for the feasibility of detecting primordial non-Gaussian 
signatures could be expressed as {Cs{0)) > {C^{0))Gauss- Unfortunately as we will see for 
the models we work with here this is far from being the result primordial features 

cannot emerge. It is worth mentioning that recent analyses of the three-point function and 
the skewness from COBE data |28| are also consistent with quasi-Gaussian fluctuations. 
Now let us turn to the examples. 

3 Worked Examples 
3.1 Natural Inflation 



To begin with let us consider the Natural inflationary scenario. First introduced in ^ 
this model borrows speculative ideas from axion particle physics Here the existence 

of disparate mass scales leads to the explanation of why it is physically attainable to have 



potentials with a height many orders of magnitude below its width [^], as required for 



successful inflation where usually self coupling constants are fine-tuned to very small values. 

This model considers a Nambu-Goldstone (N-G) boson, as arising from a spontaneous 
symmetry breakdown of a global symmetry at energy scale / ~ mp, playing the role of 
the inflaton. Assuming there is an additional explicit symmetry breaking phase at mass 
scale A ~ tugut these particles become pseudo N-B bosons and a periodic potential due 
to instanton effects arises. The simple potential (for temperatures T < A) is of the form 
r(0)=A4[l + cos(0//)]. 

For us this constitutes an axion-like model with the scales / and A as free parame- 
ters. It is convenient to split the parameter space into two regions. In the f » mp zone the 
whole inflationary period happens in the neighbourhood of the minimum of the potential, as 
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may be clearly seen from the slow-rolling equation |^ | cos(0//)| < 247r/(247r + {nip/ fY) 
which is only violated near 0end — tt/, and from the small value of the steepness |X| = 
tan(^) for (p smaller than (pend Thus by expanding V{(j)) around the minimum it is 
easy to see the equivalence between this potential and the quadratic one V ~ m^(0 — 7r/)^/2 
with w? = h!^/ p. We have already studied the latter in Ref. [|To[] and we will add nothing 
else here. On the other hand, let us consider the other regime, where / < nip. Reheating 
temperature considerations place a lower limit on the width / of the potential. For typical 
values of the model parameters involved, a temperature Trh ^ lO^GeV is attained. GUT 
baryogenesis via the usual out-of-equilibrium decay of X-bosons necessitates instead roughly 
Trh ~ lO^^'GeV (the mass of the gauge bosons) for successful reheating |Q. Thus the final 
temperature is not high enough to create them from the thermal bath. Baryon-violating 
decays of the field and its products could be an alternative to generate the observed asym- 
metry if taking place at Tpn > lOOGeV, the electroweak scale. This yields the constraint 
/ > O.Snip [0], implying n > 0.6 (see below). Attractive features of this model include the 
possibility of having a density fluctuation spectrum with extra power on large scales. Actu- 
ally for / < 0.75 nip the spectral index may be accurately expressed as n ~ 1 — mp/87r/^ 
|35| . This tilt in the spectrum as well as the negligible gravitational wave mode contribution 



to the CMB anisotropy might lead to important implications for large scale structure. [p6 



Let us take now / ~ 0.446 mp corresponding to n ^ 0.8. Slow-rolling requirements are 
satisfied provided the accelerated expansion ends by (pend — 2.78/, very near the minimum of 
the potential. Furthermore, the slow-rolling solution of the field equations yields the value 
of the scalar field 60 e-foldings before the end of inflation sin(06o/2/) — exp(— 15mp/47r/^) 
where we approximated (j)end — tt/. 

We find X|o - 4mpX^o = {nip/jy[2 + sm\<P,o/2f)]{l - sm\<P,o/2f))-' ~ 2(mp//)2. 

Two years of data by COBE ||13|] are not yet enough to separately constrain the amplitude 
of the quadrupole and the spectral index. A maximum likelihood analysis yields Qrms-ps = 
17.6 exp[0. 58(1 — n)\^K. By making use of Eq.(|^) for the rms quadrupole, COBE results 
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constrain the value of the free parameter A ~ 1.41 x 10 ^mp. We find S ^ 3.9 x 10 ^, a 
rather small signal for the non-Gaussian amplitude of the fluctuations. 

3.2 Intermediate Inflation 

We will now study a class of universe models where the scale factor increases at a rate 
intermediate between power-law inflation -as produced by a scalar field with exponential 
potential - and the standard de Sitter inflation. In Ref. p8| Barrow shows that it is 
possible to parameterize these solutions by an equation of state with pressure p and energy 
density p related by p + p = 7p^, with 7 and A constants. The standard perfect fluid 
relation is recovered for A = 1 leading to the a{t) ~ e^^"f^ {Hinj constant during inflation) 
solution of the dynamical equations when the spatial curvature A; = and 7 = 0, while 
a{t) ~ t^/^'^ for < 7 < 2/3. This non-linear equation of state (and consequently the two 
limiting accelerated expansion behaviours) can be derived from a scalar fleld with potential 
V = Vo exp(— ^370). On the other hand for A > 1 we have a{t) oc exp{At-^) (intermediate 
inflation) with A > and < / = 2(1 — A)/(l — 2A) < 1 and again in this last case it is 
possible to mimic the matter source with that produced by a scalar fleld 0, this time with 
potential 

with /3 = 4(/-i - 1). 

The equations of motion for the fleld in a = Friedmann universe may be expressed 
by 3H^ = 87T{V{(f)) + 0^/2)/mp and (p + 3if(/) = —V. Exact solutions for these equations 
with the potential of Eq.(0) are of the form |3|, |9l: H{(j)) = A/(A/?/47r)^/^(0/mp)-^/2 g^^^ 
0(t) = (A/3/47r)^/2t-^/^mp. Solutions are found for all > but only for 0^ > (/?^/167r)m|, 
we get a > (i.e. inflation). In addition f3 > 1 is required to ensure that the accelerated 
expansion occurs while the scalar fleld rolls (not necessarily slowly) down the potential, in 
the region to the right of the maximum (as it is generally the case). 
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From the full potential (|^ we may compute the value of the dimensionless skewness S 
(for convenience we will be taking the field normalized in Planck mass units from here on) 



S = 0.17Q 



f5{(3 - 4) ^ 4(3\(3 - 1) - 192n(3\(3 - 6)0' 



(8) 



02(487r02_/?2)2 

evaluated at ^ 06o- In Eq.(|]) we took I{n) ^ 4.5 (from Eq.(|^)), that is the case for the 
specific examples we discuss below. A plot of this quantity as a function of = 0^— /J^/iSvr > 
(the value of the squared of the field beyond the minimum allowed) for different values of 
f3 is given in Fig. 1. Note that both positive and negative values of S are therefore allowed 
just by modifying the choice of j3. Another generic feature is the rapid decrease of the non- 
Gaussian amplitude for increasing values of the field beyond P/^/16n (i.e. d > 0). Clearly 
this is because for large we approach the slow-roll region where the steepness becomes 
increasingly small. 

Similar calculations may be done for the spectral index. The explicit expression of n 
calculated from our potential (^ is complicated and not very illuminating. Fig. 2 illustrates 
the variation of n as a function of d (i.e. the scale dependence of the spectral index) for 
different values of the parameter f3. 

These two figures show that for acceptable values of the spectral index very small ampli- 
tudes for S are generally predicted. As an example we consider (3 = 1.2. This ansatz yields 
a negative S being 06o — 0.37 the value of the field that maximizes \S\ ~ 2.25Q. We show 
in Fig. 3 the form of the infiaton potential (|^ for this particular (3. We see that for the 
scales that exit the Hubble radius 60 e-foldings before the end of inflation |Q the value of 
the infiaton, 06o, is located in the steep region beyond the maximum of the potential. The 
value of the spectral index associated with this choice of /3 is n ~ 1.29 (a 3% below the 
first -order result n ~ 1.34). This value in excess of unity for the scales under consideration 
yields a spectrum with less power on large scales (compared with a Harrison-Zel'dovich one) 
making the long wave length gravitational wave contribution to the estimated quadrupole 
subdominant. The slow-roll parameters for this scale are e = 0.16, 1] = 0.37 and C, = 1-01 
n| . While the first two are smaller than one, ^ is not and so cannot be considered an 
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expansion variable on the same footing as e and 77. Terms proportional to ^ (and therefore 
non negligible) in Eq.(^ are those not included in the second-order analysis done for the 
first time in Ref.||23|. Taking the C05i? normalization we finally get S ^ —4.4 x 10~^. 

Let us now consider (3 = 7. Now the potential falls to zero as a power-law much more 
rapidly than in the previous case. If we take n ~ 0.8 we see from Fig. 2 that d ^ 6.05 
(0go — 2.65) corresponding to the slow-rolling region of the potential [Q. For the parameters 
we find the following values: e = 0.13, 77 = 0.17 and ^ = 0.27. From these we see that the 
first-order result (n ~ 0.86) is 7 % above the full second-order one. In this case we get 
S ~ 0.50Q. Now gravitational waves contribute substantially to the detected quadrupole. 
Actually we have ~ 1.99 |3|, |9|, |2| (while we would have had ~ 1.89 up to first 

order). Thus the estimated quadrupole should be multiplied by a factor (1 + 1.99)^^/^ to 
correctly account for the tensor mode contribution. Finally we get S ~ 7.5 x 10^^. 

3.3 Polynomial Potential 

We are interested in considering a potential of the form 

V^(0)=^(^0'+f0' + f0')+K) (9) 

where for convenience is written in Flack mass units and A, a and (3 are dimensionless 
parameters. Translation invariance allows us to omit the linear (p contribution to V. A 
detailed analysis of a potential of the form (H) was done by Hodges et al. 0. Parameter 
space diagrams were constructed and regions where non-scale invariance was expected were 
isolated. Here we will just summarize what is necessary for our study. 

Taking j3 > 8a'^/9 ensures that = is the global minimum and therefore Vq = 0. If 
we further require (3 > a"^ then no false vacua are present. Scalar curvature perturbations 
are conveniently expressed in terms of the gauge-invariant variable ( [Q, The power 
spectrum associated with it, assuming slow-roll evolution of the scalar field in the relevant 
region of the potential, is given by P^^'^{k) oc evaluated at horizon crossing time. 

Equivalently P^^'^{k) oc V^^"^ / (mpV) which suggests that regions of the parameter space 
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where the slope of the potential goes through a minimum or a maximum will be of interest 
as far as broken scale invariance is concerned. The presence of these extrema in V is 
guaranteed by taking a"^ < (3 < 4a;^/3, with a < in order for the scalar field to roll 
down the potential from the right. We will thus concentrate our study in the vicinity of 
an inflection point (pj, i.e. near the curve /3 = (but still (3 > a^). In this limit we have 
0j ~ — ^ and the number of e-foldings taking place in the region of approximately constant 
slope about 0/ is given by = —na^/ \J'i{l3 — a^) 0. Fixing = 60 as the number required 
for achieving sufficient inflation, the interesting parameters ought to lie close to the curve 
^ ~ + 7r2aV(3A^^) ~ + ya^ with ?/ ^ 9.1 x 10'^ 

We find X|q — 4XgQ = 192 y (9?/ — a~'^)/[a^(6?/ + a"^)^] where, again, we are normalizing 



the field in Planck mass units, y ~ 9.1 x 10 and we have evaluated the field for 



^60 



a 
' 2 • 



Variation of a in the allowed range results in both positive and negative values for S. 
Clearly (3 > lOa^/9 yields S > 0. Although in this region of the parameter space the 
value of a that makes S maximal corresponds to a = —5.44, this value confiicts with the 
requirement (3 < 4q;^/3 for the existence of an infiection point. We take instead a = —3.69 
which makes S the largest possible one and at the same time agrees within a few percent with 
approximating (pf ~ — |. Then /3 ~ + ya^ ^ 15.92. This guarantees we are effectively 
exploring the neighbourhood of the curve (3 = a"^ in parameter space. A plot of the potential 
for these particular values is given in Fig. 4. The slow-roll parameters in this case have 
values: e = 5.94 x 10-^ r] = 5.84 x 10"^ and ^ = 1.32 x 10"^ By using Eq.(|) we get 
n ^ 0.99. 

Let us consider now the case where a"^ < (3 < 10q;^/9. This choice of potential parameters 
leads to iS < 0. The non-Gaussian signal gets maximized for a = —2.24 and (3 = 5.13. Fig. 
5 shows the potential in this case. Note the resemblance between this form of the potential 
and that of the hybrid |^5| model Vq + m^cfP' /2 for Vq dominating and in the vicinity of the 



flat region. In that case the same sign of the dimensionless skewness ||T0| and blue spectra 
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n > 1, were predicted. For this value of a we get: e = 9.33 x 10 = 6.76 x 10 ^ and 
^ = 2.76. We get n ~ 1.01. 

We can now use Eq.(^) for the quadrupole to find the overall normalization constant A 
of the potential. Bennett et al. [|T3| get a best fit Qrms-ps = (17.6 ± 1.5)fiK for n fixed to 
one, as it is our present case to a very good approximation. Thus, for a = —3.69 we get 
A = 2.87 X 10-^2 and 5 ~ 1.1 x 10"^ while for a = -2.24 we get A = 5.88 x lO'^^ g^^^ 
S ~ -1.9 X 10-6. 

The above two values for n show that the departure from scale invariance is actually very 
small (in fact, this is because we are exploring a very narrow range of scales). Note also the 
relatively large value of ^ in the last case compared with e and 1]. This tells us that the terms 
proportional to ^ are non negligible in general. Also the rather small amplitudes for S agree 
with previous numerical analyses ^ in which adherence to the correct level of anisotropics 
in the CMB radiation under the simplest assumptions of inflation, like slow-rolling down 
with potential (^, practically precludes any observable non-Gaussian signal. 

4 Conclusions 

In the present paper we have presented explicit calculations (in the frame of some well mo- 
tivated inflationary models) of the dimensionless skewness S predicted for the large angular 
scale temperature anisotropics in the CMB radiation as well as the evaluation of the pri- 
mordial spectral index of the density perturbations originating these anisotropics. These 
computations were performed to full second order in perturbation theory. In all three mod- 
els (even in the case of the polynomial infiaton potential where more parameters were at our 
disposal) very low values for the non-Gaussian signal were obtained. 

In fact, the explicit values for S were found generically much smaller than the dimen- 
sionless rms skewness calculated from an underlying Gaussian density field and are therefore 
hidden by this theoretical noise, making experimental detection impossible. 

One may try to resort to many-field models in the hope of shifting the non-zero S window 
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to larger values. In this respect a potential of the form V{(j)i) ~ exp(X]j Aj0j) is likely to 
do well the job 0. In this case the resolution of a set of coupled Langevin-type equations 
for the coarse-grained fields (suitably smoothed over a scale larger than the Hubble radius) 
should be faced. In contrast to previous many-field analyses where one of the fields was 
assumed to dominate at a certain stage, in our case (for the aim of computing the three- 
point function) we need to make a second-order perturbative expansion in 6(f)i around the 
classical solutions but keeping V{(j)i) fully dependent on all the fields. Non-Gaussian 

fluctuations can indeed be generated within this model and thus the prospect of getting a 
non-negligible value for the dimensionless skewness should be tested. This is the subject of 
our current research. 
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Figure Captions 

Fig. 1 : Dimensionless skewness S (in units of the quadrupole Q, see Eq.(^) as a function 
of d = 0^ — jS'^/lGiT, the value of the squared of the field beyond the minimum allowed, 
for different choices of f3. Curves from bottom to top correspond to (3 from 1 to 8 (left 
panel) . 

Fig. 2 : Spectral index (as calculated from Eq.(|^)) as a function of d. Now curves from top 
to bottom are those corresponding to /? from 1 to 8 (right panel). 

Fig. 3 : Infiaton potential V{(j)) for the parameter choice /3 = 1.2 (/ = 0.77). The field is 
taken in units of mp, while the potential is normalized in units of (lO^^nipA^^^y. 

Fig. 4 : Polynomial potential V{(f)) as a function of the scalar field (0 is taken in Planck 
mass units). The overall normalization parameter A is taken to be one; for a = —3.69 
and (3 = 15.92 (left panel); 

Fig. 5 : Polynomial potential V{(j)) as a function of the scalar field {(j) is taken in Planck 
mass units). The overall normalization parameter A is taken to be one; for a = —2.24 
and P = 5.13 (right panel). 
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